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NONLINEAR STOCHASTIC WAVE EQUATIONS: BLOW-UP OF 
SECOND MOMENTS IN L 2 -NORM 



The paper is concerned with the problem of explosive solutions for 
a class of nonlinear stochastic wave equations in a domain T> C R d for 
d < 3. Under appropriate conditions on the initial data, the nonlinear 
term and the noise intensity is proved in Theorem 3.1 that the L 2 - 
norm of the solution will blow up at a finite time in the mean-square 
sense. An example is given to show an application of the theorem. 



1. Introduction. Consider the Cauchy problem for a nonlinear wave equa- 
tion 



where b\ = J^, V 2 is the Laplacian operator, and the functions f,g and h 
are given such that the Cauchy problem (1.1) has a unique local solution. 
It was first shown by Keller [6] in 1957 that, for a certain class of nonlinear 
functions f(u), the solution of equation (1.1) becomes infinite or explodes at 
a finite time, provided that the initial data satisfies appropriate conditions. 
His result was later generalized by Glassey [4] and others. Since then, it has 
become known that solutions to more general nonlinear hyperbolic equations 
may develop singularities in finite time [5], physically manifested as shock 
waves, tsunami or explosion. It is therefore of interest to examine the effect 
of a random perturbation to equation (1.1) on the existence of an explosive 
solution. To study this type of problem, it is necessary to employ some 
analytical tools from the theory of stochastic partial differential equations 
(see, e.g., [3]). 
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In our previous papers [1, 2], we considered the local and global solutions 
of a certain class of nonlinear stochastic wave equations and their asymptotic 
behavior. As an example of the nonexistence of a global solution, we showed 
that, under some explicit conditions on the initial data and the noise term, 
the solution of a cubically nonlinear wave equation perturbed by an additive 
noise will blow up at a finite time in the mean-square sense. In the present 
paper, it will be shown that this blowup result can be extended to a general 
class of nonlinear stochastic wave equations given by (2.2). Related to our 
work, we should mention an earlier paper of Mueller [7] who investigated 
the global solution of the stochastic wave equation 



for d = 1,2, where W(x,t) is a Wiener random field. He found a growth 
condition on a(u), so equation (1.2) has a unique global solution. But the 
existence of explosive solutions was not considered there. Notice that, in 
contrast with our equation (2.2), the equation (1.2) is a random perturba- 
tion of a linear wave equation. However, its noise amplitude a(u) may be 
unbounded. 

In what follows, we shall first recall some basic results for nonlinear 
stochastic wave equations in Section 2. Then, in Section 3, we prove the 
main theorem (Theorem 3.1) which states that, under some sufficient con- 
ditions on the initial data and the noise term, the solutions to a class of 
nonlinear stochastic wave equations whose L 2 -norms will explode at a finite 
time in the mean-square sense. Finally, in Section 4, we apply the theorem to 
a two-dimensional problem to obtain a set of explicit conditions for explosive 
solutions. In passing, it should be pointed out that we have only obtained 
the blowup result in a mean-square-L 2 sense. The question concerning the 
existence of explosive solutions, in the almost sure sense and at each point 
in space, is a challenging problem which remains open. 

2. Preliminaries. Let T> C M. d be a bounded domain with a smooth bound- 
ary &D. We set H = L 2 (T>) with the inner product and norm are denoted 
by (•,•) and || • ||, respectively. Let H 1 = H 1 ^) be the L 2 -Sobolev space of 
first order and denote by Hq the closure in H 1 the space of C 1 -functions 
with compact support in D. Let W(x,t) be a continuous Wiener random 
field defined in a complete probability space (Q,J-,P) with a filtration J~t 
([3], p. 38). It has mean EW(x,t) = and covariance function defined by 
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EW(x, t)W(y, s) = (t A s)r(x, y) 



x,yeV, 



where (tAs) = min(i, s) for < t, s < T. In this paper, the spatial correlation 
function r(x,y) is assumed to be bounded and continuous for x,y £ D and, 
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if D is unbounded, 

(2.1) / r(x,x) dx < oo. 



(2.2) 



Let /(//) and cr(/i,£,x,i) be two nonlinear continuous functions for /igK, 
£ € K d , x (zT> and t > 0. We consider the initial-boundary value problem for 
the nonlinear stochastic wave equation: 

' d?u{x, t) = (c 2 V 2 - a)u + f{u) 

+ a(u,Du,x,t)d t W(x,t), t>0, 
u(x,0)=g(x), dtu(x, 0) = h(x), x£V, 

I u(;t)\ & p = 0, 

with g G H 1 ,h G if, where c and a are positive parameters, and D = d x 
denotes the gradient operator. 

To consider (2.2) as an Ito equation in H, we set ut = u(-,t),vt = v(-,t) 
and so on, and rewrite it as 

{du t = v t dt, 
dv t = [(V 2 - a)u t + f(u t )] dt + dM t (u), < t < T, 
uo = 9, vo = h, 
which can be written as the integral equation 

{u t = u + v s ds, 
J % ft 
v t = v + ^ (V 2 -a)u s ds + f(u s )ds + M t (u), 

where we set 

(2.5) M t {u)= [ a s {u s )dW s , 

Jo 

and (Tt(u) = cr(u, Du, - ,t). The stochastic integral (2.5) exists under the above 
conditions on r(x,y) (see [1], the Appendix). 

Suppose that the nonlinear functions /(/i) and a(fj,,£,x,t) are locally 
Lipschitz continuous with a polynomial growth. In addition, assume that 
d < 3 and there exists an energy bound for equation (2.2): 



(2.6) 



f (v s ,f(u s ))ds+ [ TrQ s (u s )ds 
o Jo 

< c\e(ut; v t ) + C2 / e(u s ; v s ) ds + C3 for < t < T, 
Jo 

where c±, C2, C3 are some constants with c± < 1; e(u; v) is the energy function 
defined by 

(2.7) e(u;u)=c 2 ||£>u|| 2 + a||u|| 2 + |M| 2 for u £ H 1 , v £ H, 
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and 

(2.8) TrQ t (u) = / r(x,x)a 2 (u,Du,x,t)dx. 

Jv 

Then, by [2], Theorem 3.1, it has a unique continuous global solution ut E Hq 
with dtUt E H, t G [0, T], for any T > 0. Moreover, the following energy 
equation holds 

e(u t ;v t ) =e(u ;v ) +2 / (v s , f(u s )) ds + / TrQ s (u s )ds 

(2.9) ( 70 j ° 

+ 2 / (u s , <r s (u s ) dW fl ) a.s., for < t < T. 
Jo 

We remark that, in order to allow a polynomial nonlinearity, the restriction 
d < 3 is needed to prove the theorem cited above. Now, we define 

(2.10) ct>(t) = \E[{u tl u t )} = \E\\u t \\ 2 . 

Then, similar to [1], Lemma 3.1, we can compute its first two derivatives 

4>'{t) = E[(u t ,v t )} = (u ,v ) + E AlKf " c 2 \\Du s \\ 2 

(2.11) J ° 

- a\\u s \\ 2 + (u s ,f(u s ))}ds, 

(2.12) cp"{t) = E{\\v t \\ 2 - c 2 \\Du t \\ 2 - a\\u t \\ 2 + (u t , f(u t ))}. 

On the other hand, without the energy bound, one can only assert the 
existence of a unique local solution. That is, there exist a unique solution 
ut as before for t < r and a sequence of stopping times {r„} such that 
i~n T T < oo and the following holds for any n: 

rtAr„ 

UtAr n =UQ+ I V sATn ds, 

(2.13) 



•tAr n 

vtAr n = v + I (V 2 - a)u shTn ds 



tAr n 

f(u sATn )ds + M t A Tn (u). 





In particular, the limiting stopping time r may be a positive constant. It is 
clear that if P{t = oo} = 1, then the local solution becomes a global solution. 

3. Explosive solutions. Consider the local solution of the problem (2.2) 
in T> £ M. d with d < 3. We shall show that, under a set of sufficient conditions, 
the second- moment i?||ut|| 2 of the local solution will explode at a finite time 
in the L 2 -norm. To this end, we first impose the following conditions: 

(Al) The function /(/i), for \x G R, is locally Lipschitz continuous and may 
be of polynomial growth. 
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(A2) For /igR, £ G R d , x G D and £ > 0, the function cr(/x, £, x, t) is contin- 
uous in x,t and it is locally Lipschitz continuous in \x and £. 

(A3) / °° f v r{x, x)q(x, t) dx dt < oo, where q(:r, t) = sup^ ^ cr 2 (/x, £, x, i), for 
/i G R and £ G R d . 

We will state and prove the following theorem concerning explosive solu- 
tions. 



Theorem 3.1. Let conditions (Al)-(A3) hold true. Suppose that, for 
uq G H , Vo G H , the problem (2.2), or the system (2.3) has a unique contin- 
uous local solution ut G Hq with dtUt G H . Then there exists a time T e > 
such that 

(3.1) lim E\\u t f = +oo, 

provided that the following conditions are satisfied: 
(Bl) The inner product (uo,vq) > 0. 

(B2) (F(no),l)>^|e(n ;t;o) + ^ J^r(x,x)q(x,t) dxtftj 

(B3) (u,f(u)) > ^(F(u),l), for any bounded continuous function u on T>, 
where F(u) = Jq f((i)d/j, and (F(u),l) = J v F(u(x)) dx. 

Proof. We shall prove the theorem by contradiction. Suppose the so- 
lution is global so that, for any T > 0, 

(3.2) E\\u t \\ 2 <oo forO<t<T. 

Define (f>(t) = ^E\\ut\\ 2 by (2.10). We will show that 4>(t) becomes unbounded 
in finite time, or the inequality (3.2) is false. Alternatively, for any A > 0, it 
suffices to show that 

(3.3) if;(t)=<t>- X (t) 

goes to zero in finite time. To proceed we first compute the first two deriva- 
tives of ip to obtain 

(3.4) i>\t) = -\^ x+i \t)^'(t) = -\r ix+1) mi(u t ,v t )] 

and 

W{t)f 



r(t) = -\r {x+1 \t) </,"(*)- (a +i) 



(3.5) 



< -A0-C A+1 )(t){^'(t) - 2(A + l)E\\v t \\ 2 }, 
where use was made of (2.11) and the inequality: [</>'(i)] 2 < -E'||ut|| 2 £ , ||i>t | 



P.-L. CHOW 

Notice that 

r \v s J(u s ))ds= f [F(u t ) - F(u )]dx = (F{u t ) - F(u ),l). 
Jv 

From this equation together with the energy equations (2.9) and (2.7), we 
can deduce that 

E\\v t || 2 = (e(«o; v ) - 2(F(u ), 1) + E f TrQ s {u s ) ds\ 

(3-6) 1 J ° 5 

+ E{2(F(u t ),l) - c 2 \\Du\\ 2 - a\\u\\ 2 } for < t < T. 

Upon substituting equations (2.12) and (3.6) into the inequality (3.5) and 
after some simplification, we can obtain 



<(t) <-A^)- (A+1) (t)( 



2A 



2(F(u ), 1) - e(u ;v ) - E f TrQ s (u s ) ds 

Jo 



(3.7) + (2A - l)E(c 2 \\Du\\ 2 + a||u|| 2 ) 

+ E[(u t ,f(u t ))-4X(F(u t ),l)]\ 
By choosing A = 1/2 in (3.7) and making use of condition (A3), we get 
i'"(t) < -I^ 3 / 2 (t)| 2(F(u ),l)-e(u ;v ) 

(3.8) — / / r(x,x)q(x,t) dx dt 

Jo Jv 

+ E[(utJ(ut))-2(F(ut),l)]}, 

which, in view of conditions (B2) and (B3), implies that 

(3.9) ifj"(t)>0 forO<t<T. 

In the meantime, by setting A = 1/2 in (3.3) and (3.4) and taking (2.10) and 
(2.11) into account, we can obtain 

(3.10) ^(0) = ^>0 

ll n o|| 

and 

/2 

(3.11) ^(o) = _-JL_(u 0>wo )<0, 

I poll 

by invoking condition (Bl). Take T > T = -^(0)/V>'(0) = \\u \\ 2 /{u ,v ). 
We can easily deduce from (3.9)-(3.11) that the function ip(t) > is convex, 
strictly decreasing and it approaches zero at a time T e <Tq. The theorem is 
thus proved. □ 
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Remarks. 

(1) In the theorem, the domain T> may be unbounded or M. d , and the Lapla- 
cian V 2 in equation (2.2) may be replaced by an uniformly elliptic self- 
adjoint operator J2ij=i d Xi (aij(x) d Xj ) with smooth coefficients = a^, 
and 

(2) the noise term a(u, Du, x, t) dtW(x, t) can be replaced by Ya=i a i{ u > Du, x, 
t) dtWi(x, t), where ai{u, Du, x, t) are similar to a(u, Du, x, t), and Wi(x,t) 
are Wiener random fields with covariance functions EWi(x, t)Wj{y, s) = 
(t As)rij(x,y), for i,j = 1, . . . ,n. In this case the expression r(x,x)q(x,t) 
in condition (B2) should be changed to J27=i r u( x , x)qi(x, t), where qi{x, t) - 

4. An example. As an example, let us consider the following problem in 
the half plane V = {x = (xi,x 2 ) G M. 2 \xi > 0}: 

f d 2 u(x, t) = (c 2 V 2 - a)u + Xu^- 1 

+ o-q tan-^l + \Du\ 2 )e~ vt d t W(x, t), t>0,x&V, 



(4.1) 



(3 1 
u(x,0) = , d t u(x,0)- 



1 + |x| z 1 + \x\ z 

[u(x,t)\ Xl=0 = 0, 

where W(x,t) is a Wiener random field with the covariance function 

(4.2) r(x, y) = tq exp{— p(x ■ y)} ioix,y£V. 

In the above equations, A, j3, p, v, <to, vq are all positive constants, p > 2 is an 
integer, and the dot product x ■ y = x\y\ + x 2 y2- 

We can apply Theorem 3.1 to show that the equation (4.1) has an ex- 
plosive solution by choosing the above constants properly. To this end, it 
is easy to check that the conditions (Al)-(A3) are satisfied. It remains to 
verify conditions (B1)-(B3). Referring to equation (4.1), we have 

dx 1 

, i I9N9 = nP f > 



(4.3) 


(u ,v ) = (3 / 
Jv 


(4.4) 


(«,/(«))= A / 
Jv 


(4.5) 


ru 

F(u)=\ / 
Jo 


and 




(4.6) 


(F(u),l) = 



u 2p dx, 

d.t = 

2p 



- [ u 2 ?dx. 



In view of (4.3), (4.4) and (4.6), the conditions (Bl) and (B3) are clearly 
met. Finally, to enforce condition (B2), we must compute several integrals. 
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From (4.6), we get 

(4.7) <*■(«.).!> -££( r 



\ 2 p , Avr/^ 
ax 



+ |x| 2 / 4p(2p-l) 
Since a = o'otan~ 1 (l + \Du\ 2 )e~ ut , we have 

(4.8) g(rE,t)=a 2 S up[t a n- 1 (l + |e| 2 )] 2 e- 2 ^= (^f} 

so that 

2 



2 



(4.9) 



/°° / r(x,x)g(M)dxdi = r (^V) [°° I e- p ^- 2ut dxdt 
Jo Jv V 2 / Jo Jv 



r £7 2 7r 3 



To obtain the initial energy e(uo;vo), we need to evaluate the following 
integrals 

/a n «\ n 1 1 2 [ dx 7T 2 2 2 



v{l + \x\ 2 ) T 



and 



(4.11) 11^ = ^^-^*=^. 

By the definition of the energy function (2.7), we can obtain from (4.10) and 

(4.11) that 

(4.12) e(u ; v ) = ~{1 + (a + 2c 2 /3)/3 2 }. 

After substituting (4.7), (4.9) and (4.12) into the condition (B2) and sim- 
plifying the result, it gives the following inequality 

(4.13) A> ?<V_i){i + (a + 2cW + ^ 

The above inequality shows that, for a sufficiently large amplitude A of the 
nonlinear term, keeping all other constants fixed, the mean-square solution 
E||ut|| 2 of equation (4.1) will explode at a finite time. Alternatively, regard- 
ing (3 as a parameter and fixing all other constants, then the solution will 
also explode if initial displacement uq is large enough. 
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